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ABSTRACT 

Taken at face value, the observed properties of the central object in T| Car suggest a 
very massive, hot main-sequence star, only slightly evolved. If this is so, the star’s 
extraordinarily high steady rate of mass loss must dynamically perturb its outer envelope 
down to the iron convection zone, where the kinetic energy associated with turbulent 
convection can be directly fed into mass ejection. Runaway mass loss, triggered by either 
internal (pulsational, rotational) or external (tidal) forcing, would produce a secular 
oscillation of the outer envelope. In either case, the oscillation is potentially able to account 
for the observed -5 yr cycles of visual outbursts in T| Car, including the giant eruption of 
1843. 


Subject headings : convection— stars: individual (q Carinae)— stars: oscillations— stars: 
variables: other (luminous blue variables)- turbulence 



1. INTRODUCTION 


Eta Carinae is one of the brightest nonexplosive stars in the Galaxy. Often classed 
as a luminous blue variable (LBV), it displays visual outbursts, occasionally of large 
intensity, at intervals of ~5 yr (see the recent reviews in Morse, Humphreys, & Damineli 
1999; also Davidson & Humphreys 1997; Stothers and Chin 1997). Another oddity of this 
star is its large surrounding dust nebula, which was ejected in a massive eruption that began 
in (or shortly before) 1 843 and continued for about 20 yr. Despite being still enshrouded 
by dust, the hot central object is believed to have the following characteristics: luminosity, 

L= 5 x 10 6 L e if single, or possibly as low as 2.5 x 10 6 L 0 if double; effective temperature, 
T e = (1.5 - 4) x 10 4 K; and current rate of mass loss, -dM/dt = (0.3 - 3) x lO^Mg yr 1 . 
Although the star's ejecta are clearly enriched in helium and nitrogen, the measured helium 
abundance, Y = 0.3 - 0.4 (Davidson et al. 1986; Dufour et al. 1999), does not exceed by 
much the star's natal value, Y = 0.27. In view of these observations, and because the present 
rate of mass loss (believed to be that of the steady-state stellar wind) cannot be sustained for 
more than ~10 5 yr without dissipating the star's entire mass, T| Car is most likely a very 
young main-sequence star of 150-300 M 0 that is evolving quasi-homogeneously because of 
heavy mass loss (Stothers & Chin 1999). The a priori probability that r| Car is a main- 
sequence object is high, because only a few nonexplosive stars with L > 2 x 10 6 L 6 are 
known and all of them are quite similar in type (Humphreys, Davidson,& Smith 1999). 

What, then, could possibly be causing the violent instability? Although many 
theories have been proposed, none has yet gained general acceptance. 

Therefore, it seems worth while to adopt, as a working hypothesis, the notion that r\ 
Car is, in effect, a single, massive hot main-sequence star. Since at least the large eruption 
of 1843 is believed to have originated deep within the outer part of the stellar envelope, any 
hydrodynamical instabilities that are triggered in the hydrogen and helium ionization zones, 
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which for such a hot star must lie near the surface, can probably be counted out as a primary 
cause. 

Damineli (1996) and Damineli, Conti, & Lopes (1997) have suggested that the -5 
yr cycles in r| Car are caused by two colliding stellar winds at periastron in an eccentric 
double-star system with an orbital period of 5.52 yr. But the massive eruption of 1 843 
remains unexplained. Since the most recent, and best studied, outburst looks in many ways 
like a LBV (S Doradus-type) event (Sterken et al. 1999; Davidson et al. 1999; McGregor, 
Rathbome, & Humphreys 1999), all the outbursts can probably be treated as arising from 
instabilities deep in the outer envelope of a single star (which may, however, be triggered by 
periodic tidal forcing by a stellar companion). 

Steady mass loss due to the stellar wind in T| Car is so large that the star’s outer 
envelope must, from this cause alone, be constantly in a dynamic state (Stothers & Chin 
1997). Here, we examine more closely the effect of the enormous stellar wind on the 
envelope structure. A potential outburst mechanism is found to exist in a large region of 
turbulent convection that lies deep within the outer envelope. 

2. STELLAR MODELS 

Ionization of the heaviest elements-- especially the most abundant heavy metal, iron-- 
produces a large number of weak spectral absorption lines. These lines dominate the 
stellar opacity in the temperature range 10 5 - 10 6 K and furnish two local opacity peaks: a 
large peak around 2 x 10 5 K and a smaller one around 1 .5 x 10 6 K (Rogers & Iglesias 
1992; Iglesias, Rogers, & Wilson 1992). As the stellar mass increases, envelope densities 
decline, and , consequently, electron scattering becomes a larger fraction of the total opacity, 
making the two iron opacity peaks more conspicuous by contrast. Since an increase of 
stellar mass also raises the contribution of radiation pressure to the total pressure, the 
envelope becomes less stable against convection. Owing to the main iron opacity peak, 



convection manages to break out if the stellar mass exceeds -7 Above -120 the 
smaller opacity peak leads to a second, but much weaker, convection zone. These two 
convection zones have been termed "iron convection zones" (Stothers & Chin 1993 a,b). 

Although the derived threshold masses depend somewhat on chemical composition, 
rotation rate, and evolutionary age, they have been stated here for nonrotating zero-age main- 
sequence stars with (hydrogen, metals) abundances of (X, Z) = (0.70, 0.03). The radius 
fractions covered by the two iron convection zones and by the convective core are shown as 
a function of stellar mass in Figure 1. The models strongly resemble those calculated with 
the large (though erroneous) Carson (1976) opacities, which yielded a large convection zone 
in the envelope at roughly the same temperatures (Stothers 1976, Fig. 3). However, our new 
models are not very sensitive to the choice of convective mixing length, which is here taken 
to be 1 .4 times the local pressure scale height. 

To study the effects of a lower metallicity and of an evolutionary helium enrichment, 
we set also (X, Z) = (0.70, 0.02) and (0.35, 0.03). The second choice simulates a stage in 
the quasi-homogeneous evolution of a very massive star, which can be represented fairly 
well by stellar models with a homogenous composition for not too small values of X 
(Stothers & Chin 1999). Uniform rotation is treated by employing the rotational equations 
used by Sackmann & Anand (1970) and applying them to the limiting case of rotational 
breakup at the equator. 

In the case of main-sequence stars with luminosities up to 2 x 10 6 L e , mass loss 
probably has no important dynamical effect on the outer envelope. To justify this assertion, 
we use in the stellar models the standard mass-loss rates of Nieuwenhuijzen & de Jager 
(1990), which have recently been supported by the work of Lamers & Leitherer (1993) and 
Puls et al. (1996). Because, however, de Koter, Heap, & Hubeny (1997) and Crowther 
&Bohannan (1997) have derived rates that are -3 times larger for the most luminous O 
stars, we shall consider also the case where the standard rates are multiplied by a factor w = 
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3. For r| Car we adopt -dM/dt = 1 x 10' 3 M 0 yr' 1 , which has an estimated error of a factor 3 
(Davidson et al. 1986, 1995; White et al. 1994; Cox et al. 1995; Hillier 1999). 

3. DISCUSSION 

The models calculated here are fully hydrostatic, and therefore any quantities derived 
from them can only approximately represent the structure of real stars that suffer mass loss 
heavy enough to affect the outer layers. Nevertheless, we suppose that these models provide 
a sufficiently accurate zeroth-order approximation for our present purposes. The depth to 
which surface mass loss directly affects the envelope structure can be represented as a 
perturbed fraction of the total mass, 

5 M/M = | dM/dt | x dyn /M (1) 

where T dyn is the dynamical time scale of the outer envelope. Since T dyn equals simply the 
period, n, of the fundamental mode of radial pulsation, for which the nondimensional 
eigenvalue co 2 = {lltlYif R 7GM lies very close to 3 in the case of very massive main- 
sequence stars (Stothers 1992), it becomes a simple matter to compute 5M/M for any 
assigned rate of mass loss, -dM/dt, by substituting n for T dyn . Then the bottom of the 
perturbed region occurs at a mass fraction (1 - q) dyn = 8M/M below the surface, where q = 
M(r)/M. For T| Car, FI is ~1 day. 

Figures 2 and 3 show, as a function of stellar luminosity, the perturbed depth (1 - 
q) dyn for w = 1 , w = 3, and -dM/dt = 1 x 10 ~ 3 M^ yr Also plotted are the upper and lower 
boundaries of the main iron convection zone. Notice that Figure 2 includes the whole range 
of likely chemical compositions for T| Car, while Figure 3 covers all possible rates of 
uniform rotation. Although helium enrichment proves to make little difference, thus 
removing evolution as an important factor for consideration, there do exist modest increases 
of the size of the iron convection zone for greater metallicities and for faster rotation rates. 
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The surface wind from an ordinary O star with L < 2 x 10 6 L e perturbs only a tiny 
fraction of the star's mass, all of it in the radiative, outer part of the envelope. On the other 
hand, the wind's influence in r| Car extends down to 1 - q = 10 8 , which falls in the middle of 
the iron convection zone. This difference is crucial. All other types of stars that possess 
convective envelopes generate their winds in or near their photospheres, where the 
convective flux is very weak. In r| Car, however, the massive wind dissipates the outer 
layers so fast that the bottom of the wind taps directly into a highly energetic field of large- 
scale turbulence. 

Assuming the adequacy of classical mixing-length theory for estimating the gross 
properties of the turbulent convective flows (Tennekes & Lumley 1972), our present models 
predict local mean convective velocities attaining -0.4 times sound speed throughout the 
iron convection zone. This means that the convective overturning time must be of the order 
of x d Since the convective heat flux reaches 0.3 of the total flux, and the turbulent energy 
flux accounts for an additional 0.01 of the total, a large reservoir of kinetic energy is directly 
available to be fed into mass loss. It seems like a reasonable approximation to assume that 
the mass-loss rate at the stellar surface is equal to a numerical constant times the mass of 
that part of the iron convection zone that lies within the perturbed region of the star. 

Roughly then, 

-dM/dt = K( 1 -q) dyn (2) 

where K is a constant. This approximation has some justification in that (2) is then 
compatible with ( 1 ). 

Let us now suppose that as a result of some fluctuation of the turbulent flow the 
mass-loss rate begins to grow. The increased mass-loss rate will perturb more of the iron 
convection zone, thereby further enhancing the mass outflow, according to (2). As matter 
streams from the surface of the star, new material from below rises to replace it. If mass 
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loss continues long enough, the equivalent of many equilibrium outer envelopes will be 
removed. This runaway process ends either when the base of the perturbed region 
penetrates down into the underlying stable radiative layers or when the initial triggering 
disturbance dies out. Thereafter, the mass-loss rate should return to something like its 
original value. 

From Figures 2 and 3 we readily predict that, to have a modest outburst, only a small 
increase of the present rate of mass loss from T| Car is needed. On the other hand, the 
maximum possible rate consistent with our stellar models corresponds to an increase of (1 - 
q) dyn by a factor of between 10 and 10 2 , putting the base of the outer envelope well into the 
radiative layers; the rate must then have reached 10 2 to 10' 1 Mg yr \ This lies in the range 
of ri Car's probable rate during the ~20 yr period following 1843, which is variously 
estimated to have been 4 x 10 1 Mg yr -1 (van Genderen & The 1984), 2 xlO 2 Mg yr 1 
(Hyland et al. 1979), 7.5 x 10 2 M e yr 1 (Andriesse, Donn, & Viotti 1978), and -10 1 M e yr' 1 
(Davidson 1989). 

The runaway process is expected to repeat cyclically. The mean cycle time for 
arbitrarily large mass-loss rates has already been determined for massive main-sequence 
stars from quasi-static evolutionary models, and is -4 yr for the luminosity of r| Car 
(Stothers & Chin 1997). This value essentially equals the thermal relaxation time of the 
outer envelope (Tg, = Eg/L), and agrees well with the -5 yr cycle time shown by r| Car. 

Since it is also close to the 5.52 yr orbital period proposed by Damineli (1996), tidal forcing 
by a stellar companion might be both tuning and regularizing what would otherwise be a 
quasi-regular cyclicity of the type displayed by the classical LBVs. 

If our model is correct, the great eruption in 1843 must have been triggered by 
considerably larger than average turbulent fluctuations in the iron convection zone. What 
actually amplified the fluctuations is conjectural, although pulsational-mode interactions, 
angular-momentum adjustments,or even and tidal forcing may have been at work. At the 



peak of the eruption, our zeroth-order approximation for the structure of the outer envelope 
must break down near the surface. The huge mass-loss rate requires an energy supply rate 



/ GM 

1 \ 

dM 


6 L = 

— + ~ v 2 1 

— 

(3) 


V R 

2 ) 

dt 



where v ^ is the terminal velocity of the ejected matter (Forbes 1968). For T| Car we adopt v K 
= 1 x 10 3 km s' 1 (Walbom, Blanco, & Thackeray 1978; Meabum et al. 1993), M = 250 M 0 , 
R = 60 R*, and L = 5 x 10 6 L e (Stothers & Chin 1997, 1999). If the mass-loss rate actually 
attained -10 1 M 0 yr' 1 during the great eruption (Davidson 1989), then 5L must have been 
~L. This large energy drain would have caused an equivalent drop in luminosity at the 
stellar surface. At the base of the outer envelope where the mass loss begins, the luminosity 
drop would have been close to zero, because the interior layers of the star would still have 
been transporting the star’s equilibrium luminosity. 

In 1843, however, a drastic surface luminosity drop was not observed. In fact, there 
was probably a bolometric brightening by -2.5 mag (van Genderen & The 1984). 

Davidson (1989) has pointed out that the ejected material, once it was liberated from the star, 
would have radiated away its substantial internal energy. Since the total energy of the 
system (star plus ejecta) must be conserved over a short time interval, this extra radiated 
energy would have equaled, almost exactly, the energy losses sustained by the star's interior. 

Hydrodynamical modeling of so complex and extreme an object as r\ Car, involving 
detailed time-dependent interactions between mass loss and turbulence (to say nothing of 
pulsation, rotation, and tides), is not yet feasible, but the arguments made here suggest a 
model of r\ Car that may explain its instability in what is possibly the most straightforward 

way. 

This work has been supported by the NASA Climate Research Program. 
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FIGURE CAPTIONS 

Fig. 1 Spatial extent of the convective zones in homogenous hydrogen-burning stellar 
models with (X, Z) = (0.70, 0.03). Note that the radius fraction contained in the convective 
core is not monotonic with increasing stellar mass (however, the mass fraction is). 

Fig. 2— Logarithm of the mass depth below the surface as a function of the logarithm of the 
stellar luminosity. Locations of the upper and lower boundaries of the main iron convection 
zone are shown for (X, Z) = (0.70, 0.03), solid lines; (.70, 0.02), dashed lines; and (0.35, 
0.03), dotted lines. The bottom of the dynamically perturbed region is indicated in the case 
of ordinary O stars with L < 2 x 10 6 L^ (w = 1 and w = 3) and T| Car, assumed in all cases 
to be nonrotating. 

Fig. 3.— Logarithm of the mass depth below the surface as a function of the logarithm of 
the stellar luminosity. Locations of the upper and lower boundaries of the main iron 
convection zone are shown for nonrotating stars ( dashed lines) and uniformly rotating stars 
at breakup ( solid lines), with (X, Z) = (0.70, 0.02). The bottom of the dynamically 
perturbed region is indicated in the case of ordinary O stars with L < 2 x 10 6 L 0 (w = 1 and 
w = 3) and T| Car. 
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Star Light, Star Bright 

Determining the present age of the Universe depends on 
knowing how fast the distant galaxies are moving way from us 
and also on figuring out how far away they lie. As everyone 
knows, time equals distance divided by velocity. Velocities are 
measured very accurately by the Doppler red shifts of spectral 
lines in a galaxy's light. To obtain accurate distances 
astronomers have traditionally relied on a class of bright yellow 
stars known as Cepheid variables, which pulse rhythmically like 
a living heart. The slower the beat, the greater the Cepheid's true 
brightness. The distance of the Cepheid and hence the distance 
of its host galaxy are obtained by measuring the Cepheid's 
apparent brightness and inferring its true brightness from its 
measured pulsation period. In principle this is no different from 
estimating the distance of a far-off street lamp at night by noting 
how faint it is. 

Applied to astronomy, this method works well only so far as we 
can detect distant Cepheids, which is not far enough to yield an 
accurate estimate of the Universe’s age. To extend the cosmic 
yardstick ten times farther, Richard Stothers and Chao-wen Chin 
have studied the luminosities and variability of a rare class of 
stars known as luminous blue variables (LBVs), which are less 
regularly periodic but much brighter than Cepheids. A classic 
example is Eta Carinae in our own galaxy (see photo). Unlike 
Cepheids, the variability of the LBVs stems from their 
occasional violent ejections of matter into space, not from a 
pulsation. Computer models that can explain this behavior have 
been formulated by Stothers and Chin. The predicted average 
interval between major outbursts agrees with observations, and 
turns out to be shorter the brighter the LBV. Although the 
characteristic LBV periods range from four to 20 years as compared to only a few months for bright Cepheids, the long wait 
needed by astronomers to obtain galactic distances with this method is expected to ultimately provide a valuable cosmological 
dividend, the age of the Universe. 
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